Role of Bulk flow in Turbulent Convection 
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1. INTRODUCTION 

At first, we congratulate Prof. Abhijit Sen for his tremendous scientific achievements 
and for building Institute for Plasma Research. He is one of the pioneers in nonlinear 
dynamics in India. I thank him for his advise and encouragement to me throughout my 
career. We wish him and his family a healthy and prosperous future. Now we present 
some of our recent results on turbulent convection in his honour. 

Thermal convection is an important phenomenon involving instability, pattern forma- 
tion, chaos, and turbulence [1, 2]. Physics of convective turbulence differs significantly 
from that of homogenous isotropic turbulence due to the presence of the thermal plates 
and the buoyancy. Thermal plumes are generated in the boundary layer near the plate, 
after which they are transported to the bulk of the fluid. 

Scaling of large-scale quantities like Peclet and Nusselt numbers, energy spectra, etc. 
are of interest in convective turbulence. They have been modelled using the properties of 
bulk and boundary layers [3, 4, 5, 6], with the work of Grossmann and Lohse [6] (GL) 
being the most recent one. The predictions of GL theory are in close agreement with the 
results of experiments and numerical simulations of Rayleigh Benard convection (RBC). 
In this paper we derive the properties of large-scale quantities, as well as the spectra of 
energy and entropy, by applying scaling arguments to the bulk flow. We assume that 
the turbulence in the bulk is fed by the thermal plumes generated in the boundary layer. 
We will show below that significant conclusions can be drawn by focussing on the bulk 
turbulence itself. Our results based on scaling arguments are consistent with the GL 
theory, and are in reasonable agreement with earlier experimental and numerical results. 

The outline of the paper is as follows: In Sec. 2 we discuss the governing equations 
and numerical details. The scalings of large scale quantities, e.g., large-scale velocity 
and large-scale temperature fluctuations for different ranges of the control parameters 
are presented in Sec. 3. The dependence of Nusselt number and dissipation rates on Ra 



and Pr are derived in Sec. 4 and 5 respectively. In Sec. 6, we present the inertial-range 
energy and entropy spectra. In Sec. 7, we briefly describe the role of bulk motion on the 
reversals of large scale velocity in RBC. We present summary and conclusions in Sec. 8. 



2. GOVERNING EQUATIONS 

Rayleigh-Benard convection (RBC) is an idealised version of convective flow in which 
a fluid is placed between two plates. The equations for RBC under Boussinesq approxi- 
mation are 

d,u+(u-V)u = + ag6z + W 2 u : (1) 

Po 

d,0 + (u-V)0 = ^u z + kV 2 6, (2) 
V u = 0, (3) 

where, u = (u x ,u y ,u z ) is the velocity field, a and 6 are, respectively, the pressure and 
temperature fluctuations from the steady conduction state (Total temperature T = T c + 6, 
with T c as the conduction temperature profile), d is the distance between the plates, and 
A is the temperature difference. Also, z is the buoyancy direction, po is the mean density 
of fluid, g is the acceleration due to gravity, and a, v, and K are the thermal expansion 
coefficient, the kinematic viscosity, and the thermal diffusivity of fluid, respectively. 

The important nondimensional parameters of RBC are the Reynolds number Re, the 
Peclet number Pe, the Rayleigh number Ra, the Prandtl number Pr, and the Nusselt 
number Nu. They are derived based on dimensional analysis as 

u-Vu Ud 

Re = -^r~ = — , (4) 
vV 2 u v ' v 

u-ve ud 

Ra = ?g_ = ^ (6) 
vV 2 u VK 

Pr = V ~, (7) 

Total heat flux (-kVT + u z T) xy 
Conductive heat flux " {-xVT) x y 

For the Nusselt number computation, the average of the heat flux is performed over the 
horizontal planes. Under steady state, the flux over each of the horizontal planes must 
be constant. Hence, we can also compute the Nusselt number by volume averaging, i.e., 

= <-«vr+..jy» = , + M\« = , + m ^ (9) 

(-kVT)^ \ K A/ \ z / > 



where u' z = u z d/K and 9' = 0/A. It is also important to note that Pe = RePr. 



The above set of equations are nonlinear, and they do not have simple solutions. 
However, theoretical models [3, 4, 5, 6] based on scaling arguments are able to explain 
the experimental and numerical results reasonably successfully. In this paper we present 
another set of scaling arguments that relate the response parameters Pe, Re, and Nu with 
the control parameters Ra and Pr. We also compare our results with earlier experimental 
and simulation results. We point out that the scaling exponents of bulk quantities for 
both free-slip and no-slip boundary conditions are quite close. However the prefactors 
for Pe, Re, and Nu for the free-slip walls are larger compared to those for no-slip walls, 
which is due to a smaller frictional force for the free-slip walls compared to the no-slip 
walls. 

We compare our model predictions with earlier experimental and numerical results. 
For some parameters, specially for Pr = and °o, we too have performed numerical 
simulations for free-slip and no-slip walls (top and bottom ones). We assume periodic 
boundary conditions along the horizontal directions. We use "Tarang" developed by 
Verma et al. [7] for free-slip boundary conditions and "Nek5000" for no- slip boundary 
conditions. For further details of the simulations, see Mishra and Verma [8]. 

In the following sections we will describe the scaling of large-scale quantities, as well 
as that of energy and entropy spectra, as a function of Ra and Pr. 

3. SCALING OF LARGE-SCALE QUANTITIES 

In this section we describe the scaling of the large-scale velocity Ul and temperature 
fluctuation 6^. One of the generic features observed in all our numerical simulations are 
the finite amplitude of Fourier 0(0,0, In) modes [8], where the three indices indicate 
wavenumber components (k x , k y , k z ). We observe that 

8(0,0,2n)«--^-. (10) 
Inn 

The above scaling occurs due to the fact that the entropy transfer from the mode 
0(n,O,n) to the mode 0(O,O,2n) is approximately equal to entropy production due to 
the w z (n,0,n) mode, specially for lower n's (n = 1,2,3). A detailed derivation is given 
in Mishra and Verma [8]. 

The mode 6 (0, 0, In) has an important consequence on the vertical profile of temper- 
ature. The averaged temperature over horizontal planes drops sharply in the boundary 
layer, and it is approximately constant in the bulk. As shown in Fig. 1, 0(0,0,2) con- 
tributes significantly to the temperature drop near the plates. These results demonstrate 
the important role played by the 6 (0, 0, In) modes in turbulent convection. 

It is also important to note that <9(0,0,2n) do not contribute to buoyancy. The corre- 
sponding w z (0, 0, In) = because the net mass flux across any horizontal plane must be 
zero. As a result, 

u-Vu«ag0 res , (11) 
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FIGURE 1. Variation of the averaged temperature profile (solid line) along the vertical direction for Pr 
= °°, Ra = 10 8 (simulation grid size=512 3 ). The averaged temperature remains constant in the bulk, while 
it displays sharp gradient near the top and bottom plates. Dashed line represents T c + 0(0,0,2), where T c 
is the conduction profile. 



where 2 es = Q 2 — 1 9 (0, 0, 2n) | 2 (here "res" stands for residual). Using these relations 
we can write equations to estimate Ui and Q L that are 

U} U, 
ci-± = ag9 KS + c 2 v^- (12) 

U L B L A B l 

<*— = J U ^ + C4K J2 (13) 

where c\ , C2, cj, C4 are constants to be determined using the data obtained in experiments 
and numerical simulations. These constants may be weak functions of system parame- 
ters like Prandtl number, aspect ratio, etc. In the following discussion we derive scaling 
of U L and 9 L in the limiting cases: (1) Re > l;Pe > 1, (2) Re > l;Pe < 1, and (3) 
Re 1 ; Pe ^> 1 . The thickness of the boundary layers is reasonably small for Case 1 . 
However, the boundary layers cover most of the bulk for the cases 2 and 3, thus making 
the scaling of the bulk flow very relevant, specially for the limiting cases. 



Case 1: Re > l;Pe> 1 

In this regime, the nonlinear terms of both Eqs. (1,2) are much larger than the diffusive 
terms. Therefore, by matching the most dominant terms of the equations, we obtain 

A 

u-V0^-w z , (14) 
a 

u Vu w ag9 ies . (15) 



Numerical simulations [8] reveal that res ~ 9 (0, 0, 2) ~ A. Therefore, 



0l~A, (16) 



U L d 



Pe = — « VRaPr. (17) 

The proportionality constant depends on the boundary condition. Experiments and nu- 
merical simulations limit the constant to lie between 0.1 to 0.4. For this case, Eq. (17) 
implies that RaPr 3> 1 . 



Case 2: Re > l;Pe< 1 

Since Pe 1, the term (A/d)u z matches with the diffusive term in Eq. (2). The 
condition Re 3> 1 implies that in Eq. (1), the nonlinear term matches with the buoyancy 
term. Therefore, 

^u z ^KV 2 e KS1 (18) 

a 

u-Vu^age r es, (19) 

or, 

res ~ RaPr, (20) 
U L d 

Pe = — ^ RaPr. (21) 

K 

Since Pe = RePr, we obtain 

Re = Ra. (22) 
Note that 9l is dominated by 0(0,0, 2n) modes, hence 

e L « A. (23) 

The condition Pe 1 implies that RaPr <C 1 for this case. 

Case 3: Re <C l;Pe> 1 

In this case, we ignore the nonlinear term of Eq. (1), but not in the temperature 
equation, which yields 

A 

u-V0 « -u z , (24) 
a 



vV 2 u « ag 



1/2 



(25) 
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FIGURE 2. Normalized Peclet number (Pe/(PrRa)° 5 ) as a function of PrRa: simulation data for Pr=6.8 
(red box) and Pr=0.2 (blue star), and experimental data for Pr=0.7 [10]. The inset shows that the 
normalized large-scale temperature is constant with relative to PrRa. [Adopted from Verma et al. [1 1]]. 



where 



lk|0(k) 



A- 6 



1/2 



appears due to the significant contributions arising from the 



pressure gradient term [9]. Therefore, 



L ~ A. 



Numerical simulations by Pandey et al. [9] indicate that ies <C 1 (0, 0, 2) | and 



On 



Ra 5 A, 



(26) 



(27) 



a 1 1/ 2 

with 8 « 0.15. Pandey et al. [9] also show that L k |0(k)' 2 1 



0.38. As a result, 



I? 



Pe « Ra 



Ra 



0.62 



ARa ^ with C, w 
(28) 



Most of the experiments and numerical simulations of convective turbulence corre- 
spond to case 1 for which Re ^> 1 and Pe ^> 1. For this case, the scaling arguments 
are in excellent agreement with the experimental [10] and numerical simulations [11] 
as exhibited by Fig. 2. The above arguments are also in general agreement with scaling 
analysis of Grossmann and Lohse [6]. So far, no numerical or laboratory experiments 
have been performed corresponding to case 2, except for Pr = that exhibits Re = Ra 
(see Fig. 3). Regarding case 3, our scaling arguments are in good agreement with recent 
simulations by Verzicco and Camussi [12], Verzicco and Sreenivasan [13], and Pandey 
etal.[9]. 

In the next section, we will discuss Nusselt number scaling. 
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FIGURE 3. The Reynolds number Re as a function of the Rayleigh number Ra for Pr = 0. The best fit 
(dotted curve) to the data yields the scaling relationship as Re = (0.06 ± 0.23)Ra 101±0 39 . 



Using Eq. (9), 
where 



4. SCALING OF NUSSELT NUMBER 

Nu- 1 = «C) =C" e (RaPr)« 2 )y 2 (04 2 s )y 2 , 
Qe(RaPr) = 



(29) 



(30) 



Here V and t stand for the volume and temporal averages respectively. For case 1, which 
is the most relevant one for experiments, we observe through numerical simulations that 
for Pr w 1 , 



C M0 (RaPr)~(RaPr) 



-0.2 



(3D 



at least up to Ra 10 8 . Using the results of the earlier section, (u' z )l ~ (RaPr) 1 / 2 . 
Therefore, 

Nu w (RaPr) 1/2 -°- 2 « (RaPr) 03 , (32) 

which is observed in experiments [5, 14, 15, 16, 17, 18] and numerical simulations 
[12, 13, 19, 20, 21, 22] up to Ra = 10 14 or so. Note that the parameters for the above 
simulations and experiments correspond to Case 1 for which Re ^> 1 and Pe ^> 1 . 

In a recent experiment, He et al. [23] report an increase of the Nusselt number expo- 
nent to approximately 0.38 near Ra tr = 5 x 10 14 . He et al. [23] attribute the above in- 
crease to the onset of the "ultimate regime" predicted by Kraichnan [ 1 962] . Interestingly, 
the transitional Rayleigh number Ra tr corresponds to Re m Pe 0. 1 x v^R^tr ~4x 10 , 
which is close to the transitional Reynolds number for the emergence of the turbulent 
boundary layer in the flow over a flat plate, as well as for the flow past a cylinder. Ahlers 



et al. [24] reported a logarithmic profile for the temperature above Ra = Ra tr , which is in 
general agreement with the logarithmic profile for the velocity in the turbulent boundary 
layer for the flow past a flat plate. This result appears to indicate a birth of a turbulent 
boundary layer, as well as destruction of the "large-scale circulation" beyond Ra = Ra tr . 
The disappearance of the large-scale circulation beyond Ra = Ra tr has been reported by 
Ahlers et al. [24] . Chavanne et al. [25] and Roche et al. [26] report the onset of "ultimate 
regime" near Ra ps 10 12 for non-smooth plates. Their results are consistent with the fact 
that roughness of the plate can reduce the onset of turbulence in the boundary layer. 

The variation of C u q ~ (RaPr) appears to be a key ingredient for the deviation 
from Kraichnan's prediction that Nu ~ Ra 1 / 2 for the "ultimate regime". It is important 
to note that the correlation function C u q (RaPr) should flatten in the ultimate regime as 
indicated in Fig. 5. The change-over is expected to take place near Ra = Ra tr 5 x 10 14 . 

For the second case with Re 3> 1 and Pe <C 1, the large-scale Fourier modes w z (k) 
and r es(k) are in phase since u z (k) = -(Kk z d/A)6 ISS (k) [see Eq. (18)]. Therefore, 
C„ e (RaPr) « 1. Also, (u' z ) L m (0/ es ) L « RaPr. Hence 

Nu-1 « (RaPr) 2 . (33) 

Since RaPr — > 0, Nu — > 1 for this case, indicating that the convective heat transport is 
negligible, consistent with the results of Kraichnan [3]. 

For case 3 too, u' z (k) 0'(k), but overall u' z is not proportional to 6' due to the 
summation. Consequently C M e(RaPr) Ra~ 015 . Using the scaling for U L and res , we 
obtain 

Nu - 1 ^ Nu ^ C ue (RaPr)(u' z 2 )y /2 (eg) l v /2 w Ra^ 5 " 015 w Ra 032 . (34) 

Note that 8 w 0.15 [see Eq. (27)]. 

No experimental or numerical results have been reported for case 2. For case 3, 
however, our scaling arguments are in good agreement with numerical results of Pandey 
etal. [9]. 

In the next section we discuss the scaling of viscous dissipation rate e u and the entropy 
dissipation rate Eq. 



5. SCALING OF DISSIPATION RATES 



The two exact relations connecting the Nusselt number with the viscous dissipation rate 
£ u and entropy dissipation rate Eq are 

Pr 2 J 4 e u 

Nu-l = (35) 
N U = ^. (36, 
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FIGURE 4. Normalized Nusselt number (Nu/(PrRa) 27 ) as a function of PrRa. Red boxes and blue 
stars are the numerical simulation data respectively for Pr = 6.8 and Pr = 0.2 from Verma et al. [11], and 
inverted black triangles are the experimental data for Pr = 6.8 from Cioni et al. [5]. Pink circles (Kerr [27]), 
green crosses (Verzicco and Camussi [12]), black stars (Verzicco and Sreenivasan [13]), and black dots 
(Stevens et al. [21]) are the numerical data for Pr = 0.7. [Adopted from Verma et al. [11]]. 




FIGURE 5. The red circles with extended chained line depict constancy of the normalized correlation 
function C u q for the periodic box (Pr =1). Assuming that the convective turbulence becomes fully- 
developed for very large Ra, we conjecture that the normalized correlation function would become 
a constant in the ultimate regime. The second curve is the correlation function for Pr = 6.8 along 
with extended Ra~° 22 for moderately large Ra, and then a constant for the ultimate regime after some 
transitional Rayleigh number [24] . 
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FIGURE 6. Dissipation-rate correlation (C £ ) as a function of PrRa for Pr = 6.8. Both viscous-dissipation 
correlation (red boxes) and thermal-dissipation correlation (blue circles) decreases with increase of Ra. 
The best fit to the data shows the scalings C £u ~ (PrRa) -0 20 (red line) for viscous dissipation correlation 
and C £g ~ (PrRa) -0 21 (blue dotted line) for thermal dissipation correlation. 



The above equations can be rewritten as 

Nu-l = |^C £ „, (37) 
Nu = PeC £e , (38) 

where C e „ = £ u /(Ul/d) and C £g = £ e /{U L d 2 Jd). 

For case 1 (Re ^> 1 , Pe ^> 1), our numerical data indicates that C e „ ~ C £e ~ (RaPr)~ 0,2 
yielding Nu ~ (RaPr) - 3 (see Fig. 6). For case 2, the flow is Kolmogorov-like, hence 
C £u « 1 and Pe = RaPr. Therefore, Eq. (37) yields Nu - 1 « (RaPr) 2 . Using e e = 
KOl/d 2 orC £e = 1/Pe, we can deduce using Eq. (38) that Nu ~ 1 . These two deductions 
are consistent with each other. 

For case 3, Pandey et al. [9] show that C £e w Ra _a with a w 0.29, hence 

Nu = PeC £e w Ra 1_?_a w Ra 033 , (39) 

since t, w 0.38. 

The above results are consistent with the Nu scaling predicted using the large-scale 
quantities. 



6. ENERGY SPECTRA AND FLUXES 



Buoyancy acts at all length scales, hence, Kolmogorov's theory in which turbulence 
is forced at large scales may not hold for convective turbulence. We need to examine 
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FIGURE 7. Plot of the normalized kinetic spectra (left) and the entropy spectrum (right) for Pr = 6.8 
and Ra = 6.6 x 10 6 . The numerical data is not consistent with either of the two models. However, the 
upper branch of the entropy spectrum fits well with k~ 2 scaling. [Adopted from Mishra and Verma [8]]. 



the phenomenology of convective turbulence carefully. We perform our analysis for the 
three case discussed earlier. 



Case l:Pe> l;Re> 1 

When buoyancy is important, and the buoyancy term is matched with the nonlinear 
term, then 

uj/lf*ag9iuiC ue . (40) 
Assuming constant entropy cascade for 9, we obtain 

£e ~ (u L d 2 Jd)C ue ps ( Ul df/l)C ue . (41) 

Using C u q ps (RaPr)~ 0,2 , we deduce 

B k ps A(RaPr) a04 (M)- 1 / 5 , (42) 
u k ps |(RaPr) a42 (M)- 3 / 5 , (43) 



or, 



E u (k) ps d (|) (RaPr)°- 84 (M)- H / 5 , (44) 

E e (k) ps A 2 (RaPrn a4 (Mn 7/5 , (45) 

1 3 

n u (k) ps -(-) (RaPr) L26 (M)- 4 / 5 . (46) 



The prefactors are derived using the scaling of large-scale quantities derived in the earlier 
sections. The above model is called the BO phenomenology [28, 29, 30, 31]. Note that 
Eo(k) is called the entropy spectrum in literature. 



If the effective buoyancy is active only at very low wavenumbers (as in Kolmogorov's 
phenomenology for fluid turbulence), then the convective turbulence can be approx- 
imated by Kolmogorov-like phenomenology, called the Kolmogorov-Obukhov (KO) 
model. This is akin to the turbulence phenomenology for passive scalar. Here, 

E (k) « n [n M (^)]- 1 / 3 ^- 5 /3 RiA 2( Ra p r )-o.i3^-5/3 (47) 

E u (k) « [n M (/:)] 2 / 3 fc- 5 / 3 ^j(|) 2 (RaPr)°- 88 (M)- 5 / 3 , (48) 

n u (k) * iQ 3 (RaPr) 3 / 2 C eM «iQ 3 (RaPr) L3 , (49) 

According to phenomenological model proposed by Procaccia and Zeitak [30] and 
L'vov [32], BO phenomenology should be active for small wavenumbers, while KO 
phenomenology for large wavenumbers. By matching the energy fluxes of the BO and 
KO phenomenologies [Eqs. (46,49)], we can deduce the transition wavenumber called 
the Bolgiano wavenumber k%o, which is 

feo~^(RaPrn a05 *4. (50) 
a a 

Inverse of the Bolgiano wavenumber is called Bolgiano length. Thus, according to 
Procaccia and Zeitak [30] and L'vov [32], BO scaling is expected for k < k^o, while 
KO scaling is expected for k > k^o- From the above arguments, BO scaling may appear 
only for large aspect ratio boxes. 



Case 2: Pe <C l;Re > 1 

Since Aui/d ~ k0i/1 , the entropy spectrum is very steep. Consequently, buoyancy is 
active only at very small wavenumbers. Hence, Kolmogorov's theory of fluid turbulence 
is valid for this case, and we expect fc~ 5 / 3 spectrum for the velocity field. Using the 
scaling described above 

E u (k) = C K o^Q 2 Ra 2 0W)- 5 /3, (51) 
E e (k) = dA 2 (RaPr) 2 (M)~ 17/3 . (52) 



Case 3: Pe > l;Re < 1 
The momentum equation yields 

. « Vw 

ccg(e ISS )i » - T . 



(53) 




FIGURE 8. Plot of normalized kinetic spectra E"(k)k 5 / 3 (KO) and E u (k)k n l 5 (BO) for Pr = and 
Ra = 1.97 x 10 4 . The spectrum is in good agreement with the KO scaling. [Adopted from Mishra and 
Verma[8]]. 
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FIGURE 9. Kinetic spectrum E" (k) (left) and entropy spectrum E 9 (k) (right) for Pr = °o and Ra = 10 8 . 
Normalized kinetic spectrum E" (k)k 13 / 3 is flat, which is in good agreement with our predictions. Entropy 
spectrum exhibits dual branches. The upper branch corresponding to 0(O,O,2«) modes fits well with k~ 2 
curve (dotted line) [9]. 



We assume that (0 re s)i = Ra S 0i, and a constant entropy flux £q = (dlui/L)C £e = 
(Qfui/l)C Ee , which yields 

E u (k) = d (|) R^- s \kd)- l3 /\ (54) 
E e (k) = E e (k)=dA 2 (kdy l/3 . (55) 

Measurement of energy and entropy spectra using laboratory experiment is quite diffi- 
cult since the measurements are performed only at some of the points inside the container 
using thermal or velocity probes. Strictly speaking, one cannot assume Taylor's hypoth- 



esis to relate frequency spectrum measured by the probes to the wavenumber spectrum. 
In any case, the experimental results are inconclusive on the spectral exponents. 

On the numerical front, Borue and Orszag [33] reported the KO scaling for periodic 
box simulations. We believe this is so because the box size of the simulation is of the 
order of Bolgiano length. Mishra and Verma [8] performed numerical simulation for 
Pr = 6.8 and 1 for aspect ratio 2y/2, and observed inconsistency with both KO and BO 
scaling (see Fig. 7), Entropy spectrum, however, exhibits dual spectra because of the 
important role played by the (0, 0, 2n) modes [8] . Thus the spectra for Case 1 remains 
inconclusive. We need higher resolution simulations for testing this regime. 

Mishra and Verma [8] also performed simulations for zero and small Pr and found 
good agreement with the scaling for Case 2, i.e., with the KO scaling (see Fig. 8). 
Regarding Case 3, the numerical results of Pandey et al. [9] are in good agreement 
with our predictions (see Fig. 9). 

7. LARGE-SCALE MODES AND FLOW REVERSALS 

Experiments on RBC reveal that in turbulent convection, the velocity near the lateral wall 
reverses randomly. This phenomenon is called flow reversal. Researchers have attempted 
to explain flow reversals using models involving stochastic resonance, plume dynamics, 
etc. [2]. Recently Chandra and Verma [34, 35] studied this phenomenon and showed that 
the flow reversals occur due to nonlinear interactions among the large-scale modes of the 
flow. In two-dimensional geometry, during the flow reversals, the corner rolls reconnect 
and the vortices with same signs merge. The newly formed vortex has flow direction 
opposite to the original vortex. It will be interesting if the above reversal mechanism is 
at work in three-dimensional convection, as well as in dynamo. The results of Chandra 
and Verma [34, 35] demonstrate that the large-scale modes play significant role in the 
flow reversal. 

8. CONCLUSIONS 

In this paper we have presented scaling of large scale quantities, as well as the spectra 
and fluxes of velocity and temperature fields in convective turbulence. The analyses 
have been performed using the scaling of bulk quantities. The results presented in the 
paper are in good agreement with those obtained in earlier experiments and in numerical 
simulations, as well as with the predictions of GL theory. These features demonstrate 
that bulk properties are very useful in determining characteristics of convection. This 
is a significant step in modelling convective turbulence. Boundary layer near the plates 
are important in convection since the energy to the bulk turbulence is supplied at the 
boundary layer. A detailed analysis combining the boundary layer and the bulk would 
be very useful is constructing a comprehensive understanding of convective turbulence. 
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